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TWO-DIMENSIONAL NONSTATIONARY MODEL OF THE PROPAGATION 
OF AN ELECTRON BEAM IN A VACUUM 
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USSR Academy of Sciences, Keldy h Institute 
of Applied Mathematics, • ,oscow 


1. The dynamics of a relativistic electron beam In Its Inherent /4* 
electrom.agnetlc field, injected into vacuum, are considered, without 
taking account of particle collisions, within the framework of the 
Maxwollr-Vlasov equations [IJ 


5 , cUv£^-^rreh 


Here, E(t,x), H(t,x) are the electric and magnetic fields, f(t,x,p) is 

the electron pulse distribution function, p, v=Wr Is the velocity, 

/"p 1} 2 2 ^ 

W(p) = vm c +c p is the electron energy, m is its rest mass, e is the 

elementary charge, c Is the velocity of light. The electron density 

n(t,x) and current density J(t,x) are expressed by the integrals over 

Impulse space 




(3) 


The electrons are Injected from a circular plane cathode, of radius 
L along the normal to It, v;lth energy 



K = 

(4) 

and producing a current 



T - 

(5) 

the density of which 

s 


/ = <r 

(6) 


*Num.bers in the margin indicate pagination in the foreign text. 
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Is constant (y and a are assigned constants). 


There are no external fields. The plane which contains the 
cathode is assumed to be equlpotehtial and, consequently, here, the 
tangential component of the electric field equals zero. 


In cylindrical coordinates r, z, where z = 0 is the plane of the ^ 
cathode, and r = 0 is the axis of symmetry, we obtain two dimensional 
problem (l)-(3)» with 3/3<j)=0,, relative to E^(t,r,z), E^(t,r,z),H (t,r,z) 
and f= 6(P^)F(t,r,z,P^,P2) , with Initial data at t=0, 

F' ^ o , 

with the boundary conditions at z=0 




at Pe>0 

at r>L ^ p^>o 


( 8 ) 

(9) 


and with intrinsic conditions on the r=0 axis 


, pch.) = rc-pj 


( 10 ) 


The following two dimensionless parameters in the problem thus 
formulated prove to be decisive 


0 





the Injected electron energy and current. 


( 11 ) 


2. The problem formulated was solved by computer. To Integrate 
Vlasov equation (2), a macroparticle model was used. The latter con- 
sists of presentation of the electron gas as a discrete set of macro- 
particles, each of which is a group of electrons with the same coor- 
dinates and Impulse. Cf course, the electron motion is described by the 
equations , 

( 12 ) 

The electron density is calculated by the formula 


n(i x).^ 



( 13 ) 
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where summing Is carried out over all macroparticles, and S is a deltoid 
function, with a carrier on the order of the dimienslons of the calcu- 
lation cell. The current density is calculated similarly. 


The calculation region is a cylinder 0<r<r^, 0<z<z^, on the outer 
boundaries of which the following conditions are laid down 


ff. - = o 

at 

3? = ^, 

- 0 

at 



(IM 


which sim.ulates the absence of electromagnetic radiation from the outside. 


Electrons leaving the boundaries z=0, z=z^, r=r^, are excluded from 
the calculation, and those striking the r=0 axis are reflected from it, 
in accordance with (10). Quite large values of r^,z^ are chosen, so as 
not to affect the result. 



Fig. 1. Current at 
Y=2, 0=0.2, 


3. We proceed to description of the 
calculation results. First and foremost, we 
note that, in all (with respect to y» a) 
versions, practically steady state behavior 
is observed. We dwell on transient details 
below, and we now give the general character- 
istics of these behaviors. 

The tv;o schemes of movement presented in 
Figs. 1 and 2 are typical. The current lines 
for two versions, which differ in the magnitude 
of the injected current (a=0.2 and a=2, re- 
spectively), with the same energy y“ 2, are 
presented here. 

In the first case, we have a uniform, ex- 
panding beaiTi, propagating along the z axis. 

In the second case, the movement pattern is 
completely different. There is an envelope 
in the family of trajectories, the z=z^{r) 
surface, which reflects electrons ("pseudo- 
cathode''), and a nonuniform., multivelocity flow. 


The electrons flow past 
the r/z^<l cone, and a sub- 
stantial fraction of them 
returns to the cathode. 

4. A priori, it is 
clear that conditions are 
pi’Oduced in the pseudo- 
cathode behavior, for the 
development of the known two 
beam instability [2]. There 
are two counterbeams, for- 
ward and reflected, in which 
the density of the fluxes 
is high in the reversal re- 
gion and, consequently, 
high frequencies should be /II 
excited. 

On the other hand, there are stabilizing factors, steady state in- 
jection from the nearby cathode, constant current drift and a decrease 
in density with distance from the axis. 

As the calculation results show, the interaction of all these factors 
results in an oscillatory process, the amplitude of vjhlch is small, and 
the average state corresponds to the steady state behavior pointed out 
above . 

The transient effects intensify with increase in current. There- 
fore, we turn to the version with Y"2, o* 4 to describe them. S (z) 
curves on the r=0 axis at different moments of time are presented in 
Fig. 3. The points where E passes through zero correspond to the 
pseudocathode location z^. E vs. time is shewn in Fig. 4 , at a point 
located on the r=0 axis near the pseudocathode {z^0.23L), for the same 
version, Y = 2, o = 4. In this area, because of the lai'ge field gradient 
in the direction of movement of its profile, the amplitude of the os- 
cillations is com.parable to the size of the field itself. 

4 




Fig. 3* E profiles on axis at y=2, a»^. 
z 

As a result of the 
pseudocathode pulsations, 
a small fraction of the elec- 
trons periodically prove to be 
above It , and they are dumped 
In the space beyond the cathode. 

These pulsations give rise 
to transient phenomena In the 
entire flow. Electrons which 
reach the pseudocathode along 
one current line acquire various 
velocities upon reflection from 
It and, then, moving practically rectllinearly , they are located on 
curves, the shapes of which are similar to the Jet from an oscillating 
hose. The movement of such Jets basically determines the structure of 
the weak turbulent flow which arises. 

For this reason, in descriptlcn of the steady state characteristics 
of the behaviors by means of the current lines, we have represented the 



Fig. E (t) In pseudocathode area, at 
Y=2, 2 
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latter provisionally In the figures, by a dashed line which shows the 
average movement. 


It can be said that the transient nature of the process Is reduced /I 7 
to the transient nature of the pseudocathode. 



Fig. 5. Flow at y® 2, a*0.5 



Fig. 6. E and 1\. profiles on 
axis, at y= 2, o=ot5 


5. We know dwell in greater 
detail on the dependence of the aver- 
age steady state properties of the 
flow on current o • We consider ver- 
sions with injected electron energies 
Y-2. 


The results of calculation of the 
version with a*=0.5 are presented in 
Flg. 5. Although, the trajectories 
intersect in a narrow boundary layer, 
on the whole, the team remains uni- 
form, with practically constant current 
density over the cross section, which 
rises slightly at the boundary Itself. 
Compared with the o=0.2 version 
(Fig. 1), the rate of broadening of 
the beam has increased appreciably. 

Curves of the E^^Ct) and pulse 

P (z) fields along the z axis are pciver 
z 

in Fig. 6, for the same version a*0.5. 

The point where E *0, P ■=mln corre- 

z z 

spends to the minimum potential. Upon 
passing it, the flow is accelerated. 
However, electron energy W(P^) rem.ains 
less than the initial W^^ymc^. 

On the same term.s, the variant 
v’lth o=-0.65 is shown in Figs. 7 and 8. 
Uniformity is preserved only in the 
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axial region. A pseudocathode 
forms in the outer layers, ^e 
outer boundary of the injected 
beam (current line V) proves to be 
inside the flow. Eigjaty percent 
of the current moves outside it. 

The version o*0.75 is pre- 
sented in Fig. 9. The pseudo- 
cathode has formed completely and 
reached the axis. The outer layer 
of the beam becomes the inner one 
(line V). The flow ’’was turned 
Inside out". The main portion of 
the current t flows in a cone 
between lines II and V. A negli- 
gible fraction of the current (•v^JO 
leaks through the pseudocathode in 
a transient manner. A reverse 
current develops. 

As we see, a qualitative re- 
arrangement of the structure and 
shape of the flow occurs over a 
small interval of change of q 
from 0.5 to 0.75. A critical cur- ' 
rent value can be spoken of, in /I 9 

this case, a 'vO.6. 

' cr 


With further Increese in 
current, the flow pattern does not change qualitatively. The pseudo- 
cathode subsides arid becomes flatter. The reverse current fraction in- 
creases. At o«2 (Fig. 2) of the electrons return on the z«0 plane 
and, at c»4, 3S%. 

f 

In the outer zone, the flow retains the shape of a cone with the 
same angle 'v45°. Of course, the fraction of the current traveling in 
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It decreases 


6. We have reported the 
calculation results of versions 
with Injection energy y»2. At 
other Y» the quantitative char- 
acteristics of the solution are 
different » of course, but all 
the effects described are pre- 
served. Fig. 10 gives an idea 
of the presence of the pseudo- 
cathode and its height (at 
r*0), in versions with other y* 

Specific information on z« 

and 0 ^^ vs. Y also can be obtaln- 
cr 

ed, by means of analytic estimates. 

7* In the limit, as o-*«, the 
pseudocathode approaches the 
cathode without restriction, 

Z||/L^0, But, with z,ic<L, the 
finite dimensions of the cathode 
should not affect the processes 
’ in the axial region. Therefore, 
if the increase of a occurs 
through an increase in cathode 
radius L, with the density of the 
injected current preserved, for 
determination of the asymptotes 
of the solution in the axial zone, 
a unidimensional version of our 
problem, 3/3r=0, which corresponds 
to an infinite cathode, can be 
used. 



8. Simple reasoning permits on estimate cf the oritioal current 

0-„(y) tc be obtained, at which a change in type of behavior cocuro. 
cr 

We find the conditions of the existence cf a steady state, single 
velocity flow, unbounded along the z axis. 

^e motion of an electron in a steady state field with potential 
0, occurs with retention of its total entrgy, W-e^>const. Since, 

at the time of injection, all the electrons have the same energy 
Wp-Yroc , and the cathode plane is equipotential , ^l2ao”^0* 

^ ( 19 ) 

for all electrons in the entire flow. 


, I The potential satisfies the Poisson equation 

r ^ y 

' ( 20 ) 

' end it is finite at infinity. i 

f S 7 

a \ 

h- i 
6- 

jy 

|- - 5 

I ^ ■ -1 

I 5 . 

t = 

7 9 ' 

" i 


Por a 9 Ingle velooity flow, the current density 


( 21 ) 

By vlsrtue of the assumed boundedness of the beam with respect to 
r, at each s, evidently, as . This gives a basis for pre- 

senting the radial portion of the Laplaoe operator in the fom 

( 22 ) 

where a is a dimensionless, indeterminate coefficient, and R is the £21 
efl^otlve radius of the beam, such that 

(23) 


is the total current. 

With (19), (21)-(23) taken into consideration, we rewrite equation 
(20) in the fom 

On the assumption that R/z < const, it is easy to determine that a 

necessax 7 condition of the existence of a bounded solution of equation 

(24), as z -M» , is the possibility of the reduction of its right side 

to zero. Othez*wlse, either or the solution exists only in a finite 

2 

interval, since W cannot be less thu: me . 

By making the right side of (2^) equal tc zero, with fixed Iq and 
a, we obtain sn equation for V. Simple analysis shows that it has a 
real root, only in the event 

«•<*>' (25) 

This formula gives an estimate of the value of the critical cur- 
rent If it is -i good agreement with the calculation 

results (Fig. 10). In particular, this justifies the ’’simplicity*' of 
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the reanoning used 
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